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ABSTRACT 

Wc look at observational constraints on the thawing class of scalar field models pro- 
posed to explain the late time acceleration of the universe. Using the recently intro- 
duced 'Statefinder Hierarchy', we compare these thawing class of models with other 
widely studied dark energy (and modified gravity) models to check the underlying pa- 
rameter degeneracies. We put constraints on the deviations of these thawing models 
from the canonical ACDM model using a large class of observational data, e.g, the Su- 
pernova Type la data, the BAO data, the CMB data and data from the measurements 
of the Hubble parameter using red-envelope galaxies. We also forecast constraints using 
a simulated dataset for the future JDEM SNe survey. Our study shows that, although 
with current data it is difficult to distinguish different thawing models from ACDM, 
a future JDEM like mission would be able tell apart thawing models from ACDM for 
currently acceptable values of flmo- 
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1 INTRODUCTION 

One of the great mysteries of the universe today is the na- 
ture of dark energy which drives the late time cosmic ac- 
celeration. This has been confirmed by numb er of obser- 
vational results including Type-la Su pcrnovac iessetal 
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diation (CMBR) ( Komatsu et al. (2011)) as well as the 
latest surveys of the large scale structure ( Eisenstein et 
al. (2005)). There has been wide range of contrasting pro- 
posals to explain this late time cosmic acceleration. Most 
popular of these is to include a dark candidate with a large 
negative pressure (also known as dark energy) in the energy 
budget of the universe ( Bean et al. (2005), Copeland et 
al. (2006), Li et al. (2011), Padmanabhan(2003), Peebles 
& Ratra (2003), Sahni & Starobinsky (2000)). 

The simplest candidate for this dark energy is a cos- 
mological constant A. However, the well-known fine tun- 
ing and coincidence problems render A a rather unattrac- 
tive proposal from a theoretical point of view. Going be- 
yond A, scalar field models with generic features are per- 
haps the simplest alternatives to a cosmological constant. 
A large class of scalar field models including quintessence 
( Ratra & Peebles (1988), Caldwell et al. (1998), Lid- 
die & Scherrer (1999), Steinhardt et al. (1999)), tachyon 
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( Abramo & FinelU (2003), Bagla et al. (2003), Aguirre- 
gabiria & Lazkoz (2004), Copeland et al. (2005)) phantom ( 
Caldwell (2002)) and k-essence ( Armendariz et al. (2001), 
Scherrer (2004), Sen (2006), Vikman (2005)) have been 
thoroughly investigated in recent years to explain the late 
time acceleration of the universe. The advantage of these 
scalar field models is that not only can they alleviate the 
fine tuning and coincidence problem, they can also mimick 
a cosmological constant at present epoch. These models have 
equation of state as a function of redshifts. This feature is 
important to distinguish these scalar field models from cos- 
mological constant as large amount observational data from 
higher redshifts are coming in or expected to come in near 
future. 

The simplest scalar field is the one with linear potential 
having a canonical kinteic energy. In this case, the scalar 
field is frozen initially due to large Hubble damping and 
behaves like a cosmological constant with w ~ —1. As the 
universe expands and Hubble damping decreases, the field 
starts rolling and the equation of state starts deviating away 
from w = —1. But because the potential has no minima, the 
model gives a collapsing universe in future which results in 
a finite history for our universe ( Kratochvil et al. (2004)). 

More complicated scalar field models can be divided in 
to two broad classes: the fast roll and slow roll models, also 
termed as freezing/tracking and thawing models( Caldwell 
& Linder (2005)). The fast roll models have steep poten- 
tials allowing the scalar field to mimick the background mat- 
ter/radiation, and to remain subdominant for most of the 
history of the universe. Only at late times, the field becomes 
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dominant and starts behaving like a component with nega- 
tive pressure driving the acceleration of the universe. On the 
other hand, slow roll or thawing models are similar to the 
inflaton field, where the potential is sufiiciently flat. Due to 
this, kinetic energy of the field is much smaller than the po- 
tential energy and field is initially frozen due to large Hubble 
damping in the early universe. That is why initially w ~ —1, 
and the energy density of the scalar field is nearly constant 
and have a negligible contribution to the total energy bud- 
get of the universe. But as the background radiation/matter 
gets diluted due to the expansion of the universe, the back- 
ground energy density becomes comparable to the scalar 
field energy density. At this moment, the scalar field starts 
thawing and slowly evolves away from the frozen w ~ — 1 
state towards a higher value of w. However, in this case some 
degree of fine tuning is needed in order to get a suitable late 
time evolution. One of the early studies for thawing mod- 
els was by Alam et al. ( Alam et al. (2003)) where they 
studied the decaying dark energy models which have thaw- 
ing type cosmological evolution. Recently, Scherrer and Sen 
( Scherrer & Sen (2008a))have studied the thawing scalar 
field models with a nearly flat potential satisfying the slow 
roll conditions. They have shown that under such conditions, 
thawing models with variety of potentials evolve in a sim- 
ilar fashion and results in a unique analytic expression for 
the equation of state ■w{z). This result was later extended to 
the case of phantom models ( Scherrer & Sen (2008b)) and 
tachyon scalar field models ( Ali et al. (2009)). In another 
work, thawing models ( for both canonical and non canoni- 
cal scalar field) have been studied without the assumption of 
the slow roll condition ( Sen et al. (2010)). In a more recent 
work, Dutta and Scherrer ( Dutta & Scherrer (2011))have 
studied the slow roll freezing scalar field models and derived 
an analytical expression for the equation of state it; as a 
function of the energy density of the scalar field Sl^. 

There is another elegant way to explain the present cos- 
mic acceleration. This is achieved through the large scale 
modification of the gravitational sector of the theory. These 
models include DGP model ( Dvali et al. (2000)), Cardas- 
sian model ( Freese & Lewis (2002)), modified GCG model ( 
Barreiro & Sen (2004)), f(R) theories ( Starobinsky (2007), 
Hu & Sawicki (2007),Amendola et al. (2007)), galileon mod- 
els ( Nicolis et al. (2009), Chow & Khoury (2009)) as well as 
models having non local gravitational interactions ( Deser 
& Woodard (2007)). 

With current obsevational data, the concordance 
ACDM indeed appears to provide a very good fit to the 
present observations. Hence, for any given dark energy 
model (different from ACDM), it is natural to ask how to 
distinguish it from ACDM. Keeping this in mind, different 
diagnostic measures have been proposed. Om diagnostic and 
the Statefinders are two elegant examples for such diagnostic 
measures. Recently, another interesting diagnostic measure, 
^Statefinder hierarchy, Sn' (with n > 2), has been introduced 
by Arabsalmani and Sahni ( Arabsalmani & Sahni (2011)). 
It was shown that while all the 5n's remain fixed at unity 
during the entire evolution of the universe for concordance 
ACDM, they are in general time dependent for dark energy 
models. Hence they can be used not only to compare dif- 
ferent dark energy models with ACDM but also to compare 
different theoretical approaches to explain present cosmic 
acceleration. 



In the present investigation, we consider thawing dark 
energy models involving canonical scalar fields. We do not 
assume the slow roll conditions on the potential of the field. 
We first use the Statefinder diagnostic Sn to compare dif- 
ferent thawing models with other dark energy models as 
well as with modified gravity proposal such as DGP model. 
Next, we use the latest observational results involving back- 
ground cosmology to put constraints on the deviations of the 
thawing models from the ACDM behaviour. We also forecast 
constraints from future observations. 



2 THAWING SCALAR FIELD MODELS 

We assume that dark energy is decsribed by a minimally 
coupled scalar field 4> with canonical kinetic energy, governed 
by the equation of motion 



(1) 



where H is the Hubble parameter determined by the Ein- 
stein equation 



3H ^ pm + I 



(2) 



We work in the natural unit 8ttG = c = 1. Here, pm is 
the non-relativistic matter (pressure p — 0) content of the 
universe and the p0 is the energy density of the scalar field. 
We assume a spatially fiat universe, so that Qm + ^t/> = 1. 
With this, one can form an autonomous system using equa- 
tions ^ and ^ (see Scherrer & Sen (2008a)) for detail 
derivation): 



7 = -37(2 - 7) + A(2 - 7) \/3^, 
n'^ = 3(1-7)!^^(1-!^^), 
a' = -\/3A^(r- 1)^7?^^. 



(3) 
(4) 
(5) 



Here 7 = 1 -f (w^ being the equation of state for the 
scalar field), and 



A 



_}_dy_ 



(6) 
(7) 



Given the initial conditions for 7, Sl^ and A, one can solve the 
above system for any potential V{(j>) (The information about 
the form of the potential is encoded in F). As our goal is to 
study the thawing models, the scalar field is initially frozen 
due to large Hubble damping (i.e (pi ~ 0). Hence ~ — 1 
initially which gives 7 = as initial condition. We shall not 
assume slow-roll conditions for the potential V{(j})- Hence, 
the initial value Ai is a parameter in our model. This initial 
value of Ai determines the deviation of from the initial 
~ — 1 frozen state as the universe evolves. For smaller Ai, 
the deviation is negligibly small, and the scalar field behaves 
as a cosmological constant at all time. For larger values of Ai, 
100 deviates substantially from —1 as the universe evolves. 
In general, the contribution of scalar field to the total energy 
density is negligible in the early universe. Nevertheless one 
has to fine tune the initial value of in order to have its 
correct contribution at present. Hence i}^ (initial) is related 
to the ^l,jy(z = 0). With this we evolve our system from the 
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Figure 1. The panel shows the Statefinder 5*2 versus Sl^ = 

VtrnO ^^^2^ ■ Larger value of Q,m i— > 1 corresponds to the distant 
past (z >> 1) whereas, the smaller value f!m i— ^ corresponds 
to the future. fimO = 0-3 for this plot. The different solid lines 
are (Right to Left) : Thawing quintessence model with F = 0, 
r = 0.5, r = 1, r = 1.5, r = 2 respectively. The dashed lines 
corresponds to CG (1), GCG (2), DGP (3), GPL (4), to = -0.8 
(5). The region between two vertical solid lines roughly corre- 
sponds the present epoch. 
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Figure 2. Same as Fig.l but for Sz- 



intial era z = 1000 till any redshift we want. We consider 
various power-law potentials for which F is a constant. In 
recent past, Clemson and Liddle have considered thawing 
models using a multiparameter extension of the exponential 
potential ( Clemson & Liddle (2009)). 



3 STATEFINDER HIERARCHY AND 
PARAMETER PHASE SPACE 

One can Taylor expand the scale factor of the universe 
around the present era (t — to) as follows: 



a{t) = a{to) + a{to)J2 



a„{to) 



[Ho{t-to)r.^ 



where, 



in 



(8) 



(9) 



It is easy to check that —02 = g is the deceleration parame- 
ter. Similarly as is related to the Statefinder r or jerk j and 
Q4 is related to the snap s and so on. Using the equations 
(O and ((5)), one can now define the Statefinder Hierarchy 
S„ ( Arabsalmani & Sahni (2011)) as: 
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and so on. 



(10) 

(11) 

(12) 
(13) 
(14) 



It is remarkable to see that for ACDM, Sn = 1 during the 
entire history of the universe. In the rest of this section we 
use different combinations of Sn, to study the evolution of 
thawing models and compare them with other popular dark 
energy and modified gravity models. We consider the follow- 
ing models: 

• Dark Energy with a constant equation of state, w = 
Too = constant. 

• Dark Energy with a CPL parametrized equation of state 
( Chevallier & Polarski (2001), Linder (2003)), w = wq + 
Wa{l — a), where wo and Wa are constants. 

• Dark Energy with a Generalized Chaplygin Gas 
( Bento et al. (2002)) type equation of state, w = 
— ^ ^"^^ -3(1+3) 1 where As and /3 are constant. (3 — 1 
gives the standard Chaplygin gas. 

• Dvali, Gabadadze and Porrati (DGP) model ( Dvali et 
al. (2000)). 

In Fig.l, we show the behaviours of the models in the 
S2 — fim phase plane. The figure shows that in this phase 
plane which involves the second derivative of the scale fac- 
tor, the models are highly degenerate around the present 
value of fim. As we go to the higher order Sn's, one can 
remove these degeneracies among various models. This is 
revealed in the subsequent figures. Fig. 2 and Fig. 3 where 
we show the phase portrait in the ^3 — Qm and 5*4 — 
planes respectively. In Fig. 4 we show the phase portrait in 
5*4 — 5*3 plane where again the degeneracies among different 
models are broken. We should comment on 5*3 in particu- 
lar. From Fig. 2 and Fig. 4, one can see that S3 can be 
used to break the degeneracies among the thawing models 
themselves as its present day value varies substantially for 
different thawing potentials. Similarly, S4 lets us differenti- 
ate between thawing and other models. We should mention 
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Figure 3. Same as Fig.l but for 54. 
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S3 

Figure 4. The phase plane for S4VS.S3. Dashed Hnes are same 
as Fig 1. The sohd hnes (top to bottom): Thawing quintessence 
model with T = 0, T = 0.5, T = 1, T = 1.5, T = 2 respectively. 
n^O = 0.3 for this plot. The point (1,1) is ACDM. 



that although we assume QmO = 0.3 for these figures, the 
conclusions remain same for other values of SlmO- 

With future data, one can probe the background expan- 
sion at higher redshifts thereby probing the higher deriva- 
tives of the scale factor. Hence Statefinder Hierarchy will be 
quite useful to distinguish various dark energy models. We 
further study the w' — w phase plane for different thawing 
dark energy models. Here 'prime' denotes the differentiation 
with respect to log a . As discussed in the previous section, 
for smaller values of Ai , all the thawing models behave simi- 
lar to the cosmological constant. So to study the w'—w phase 
plane, we assume = 1 for which the thawing models be- 
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Figure 5. The panel shows the dynamics of scalar field in w — w' 
plane at redshifts 2 = (solid)a,nd z = I (dashed).Thc dots 
(black) represents the potentials V = 4>, (p^ , (p"^ , (p^^ (from top 
to bottom, for each z)respectively. 



have significantly different from cosmological constant. Lin- 
der and Caldwell have previously showed that the thawing 
models are constrained to lie between (l-|-«;) < u;' < 3(l-|-u;) 
region in w' — w phase plane. In Fig. 5 and Fig. 6, we show 
different thawing models in the w' — w phase plane at dif- 
ferent redshifts and compare them with the constraint pro- 
posed by Linder and Caldwell ( Caldwell & Linder (2005)) 
(LC). It can be seen that for higher redshifts, all the thaw- 
ing models are inside the LC bound. This is obvious as all 
the thawing models behave very close to cosmological con- 
stant for larger redshifts and hence satisfy the LC bound. 
But at the present epoch (z=0), for smaller values of QmO, 
thawing models with certain potentials ( as shown in Figure 
5) are outside the LC bound. But as one increases the value 
of Qmo (see Figure 6), they tend to satisfy the LC bound. 
Hence thawing models which deviate significantly from the 
ACDM behaviour at smaller redshift tend to satisfy the LC 
bound for higher values of ^ImO whereas for lower values of 
fimo, some of the potentials do not satisfy the LC bound. 



4 OBSERVATIONAL CONSTRAINTS 

In this section, we constrain the parameters in the thaw- 
ing scalar field models with the assumption of a fiat Uni- 
verse by using the latest observational data including the 
Type-la Supernovae Union2 compilation ( AmanuUah et 
al. (2010)), the BAG ( Eisenstein et al. (2005)) measure- 
ment from the SDSS ( Percival et al. (2007), Percival et 
al. (2010)), the CMBR measurement given by WMAP7 ( 
Komatsu et al. (2011)) observations, the H(z) data from 
HST key Project ( Freedman et al. (2001)) and a simulated 
dataset ( Holsclaw et al. (2010)) based on the upcoming 
JDEM SN-survey containing around 2300 Type-la Super- 
novae. 
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Type la Supernovae 

We consider the Supernovae Type la observation ( Aman- 
uUah et al. (2010) which is one of the direct probes for the 
late time acceleration. It measures the apparent brightness 
of the Supernovae as observed by us which is related to the 
luminosity distance c?l(z) defined as 



dL{z) = {l + z) 



dz' 
W) 



(15) 



With this, we construct the distance modulus jj, which 
is experimentally measured 



jj, = m — M — 5 log 



Mpc 



+ 25 



(16) 



Where m and M are the apparent and absolute magnitudes 
of the Supernovae which are logarithmic measure of flux and 
luminosity respectively. 



Baryon Acoustic Oscillation 

Another observational probe that has been widely used in 
recent times to constrain darli energy models is related to the 
data from the BAO distance measurements obtained at z — 
0.2 and z — 0.35 from the joint analysis of the 2dFGRS and 
SDSS data ( Eisenstein et al. (2005), Percival et al. (2007), 
Percival et al. (2010)). In this case,one needs to calculate 
the parameter _D„ which is related to the angular diameter 
distance as follows 



zbao 

H{ZBA0) 



dz 
H{z] 



1/3 



(17) 



For BAO measurements we calculate the ratio 
D^iz=0.35) 



D^{z=0.20) 



this ratio is a relatively model independent 



Figure 7. The Icr and 2a confidence contours in Ai - CmO plane 
for r = O.The shaded regions are constraints from SN+BAO data 
while the thick lines are constraints from SN+BAO+CMB+H(z) 
data respectively. The dashed lines are for the simulated JDEM 
data. The two vertical lines represent the WMAP7 bound on QmO 




quantity and has a value 1.736 ± 0.065. 



Figure 8. Same as Figure 7, but with F = 0.5. 



Cosmic Microwave Background 

The CMB is sensitive to the distance to the decoupling 
epoch via the locations of peaks and troughs of the acous- 
tic oscillations. We employ the "WMAP distance priors" 
given by the five-year WMAP observations. This includes 
the "acoustic scale" ^Ajthe " shift parameter " R and the 
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Potential, V = ^ 
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Figure 9. Same as Figure 7 but witii F = 1.5. 

redshift of the decoupling epoch of photons The acoustic 
scale Ia describes the distance ratio ^^y*'' 

'-<' + -'^ (18, 

where (l + z,) factor arises because Da{z*) is the proper an- 
gular diameter distance, where ra{zt) is the comoving sound 
horizon at z. 

We use the fitting function of z, proposed by Hu and 
Sugiyama ( Hu & Sugiyama (1996)) 

= 1048[1 + 0.00124(^6/1^)-°-^^*] [1 + gi{n„,h'^y'-] (19) 
_ 0.0783(n6fe^)-°'^^* _ 0.560 

1+ 39.5(f76/l2)0.763 ' f2 - ^ _^ 21.1(06/12)1-81 ^ ' 

The shift parameter R ( Komatsu et al. (2011)) is re- 
sponsible for the distance ratio ^-i^' given by 

R{z^) = s/^^oHlil + z,)Da{z*) (21) 

The constraint on the shift parameter R{zt) from the 
WMAP observations is quoted as R{z,) = 1.715 ± 0.021. 

Expansion history H(z) data 

We use new determinations of the cosmic expansion history 
from red-envelope galaxies. Stern et al. ( Stern et al. (2010)) 
have obtained a high-quality spectra with the Keck-LRIS 
spectrograph of red-envelope galaxies in 24 galaxy clusters in 
the redshift range 0.2 < z < 1.0. They complemented these 
Keck spectra with high-quality, publicly available archival 
spectra from the SPICES and VVDS surveys. With this, they 
presented 12 measurements of the Hubble parameter H(z) 
at different r edshift. The measurem ent at z = is from HST 
Key project jPreedman et al.|[200ll '). 



Potential, V = Exp(A0) 




Figure 10. Same as Figure 7 but with F = 1. 



Simulated SN data 

To look at what can be achieved in future, we also use a 
simulated dataset ( Holsclaw et al. (2010)) based on the 
upcoming JDEM SN-survey containing around 2300 SNe. 
These are distributed over a redshift range from z = to 
z = 1.7. We assume a simplified error model where the errors 
are same for all supernovae and independent of redshift. We 
assume a statistical error of cr = 0.13 magnitude, as expected 
from JDEM-like future surveys( Aldering et al. (2004)) . We 
use a ACDM model with Ho = 72 km/s/Mpc and ilmo = 
0.27 to simulate the data. 



Constraints on Model Parameters 

Thawing models are described by three parameters, namely 
Ai, r and fimo- As mentioned in the previous sections, Xi 
describes the deviation of the thawing model from ACDM 
(Ai — represents the ACDM model), and F encodes the 
information about the potential. Hence constraining these 
two parameters essentially constrain the deviation of differ- 
ent thawing models from ACDM. In Figures 7, 8, 9 and 10, 
we show the constrained region in Ai — QmO plane for differ- 
ent potentials. In doing this, we use the data coming from 
different recent observations e.g Snia, BAO, CMB as well as 
H(Z) measurements. We also use the simulated data for the 
upcoming SN-surveys like JDEM to constrain the parame- 
ter space. In these figures, we also show the WMAP7 bound 
on QmO- 

From the figures it is evident that when we consider 
the current data, there is always an upper bound on Ai. 
This upper bound is most prominent for the linear poten- 
tial. Hence this gives strong constraint on deviation from 
the ACDM behaviour for the linear potential. For Expo- 
nential or Inverse Squared potential, this deviation is not 
strongly constrained and Ai as large as 1 is allowed even at 
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Figure 11: The behaviour of S3 at present {z = 0) as a funtion of for different thawing potentials. The lines represent 
r = 0, 0.5, 1, 1.5, 2 from bottom to top. QmO = 0.24, 0.26, 0.28 (from left to right) respectively. The horizontal lines represent the bound 

on 53 at 2 = due to the JDEM simulated data. 



2 — a confidence level. Moreover if one takes into account 
the WMAP7 bound on QmO, current data always allow the 
ACDM behaviour (Ai = 0) for all the potentials. 

But this changes as one does similar analysis with 
the simulated data for future SN-surveys like JDEM ( 
Holsclaw (2010)). With this, considering the WMAP7 bound 
on ilmo, although ACDM (Ai = 0) is allowed for some val- 
ues of QmO, for other values of QmO, Ai has a lower bound 
different from zero. This clearly rules out ACDM even at 
2 — (J confidence level. For example, consider the constraints 
obtained from simulated JDEM dataset in case of linear po- 
tential (Fig. 7). The contours are shown by the dashed line. 
For n,„o ~ 0.25 (which is allowed by the WMAP7 constraint 
on Qmo), the lower bound on A; is 0.75 at 2 — cr confidence 
level. Similarly, in case of an exponential potential (Fig. 10), 
the corresponding lower bound on A; for the same QmO is 
0.85 at 2 — (T confidence level. This rules out the ACDM 
even at 2 — cr confidence level. For all the potentials that we 
considered, one can rule out ACDM at 2 — a confidence level 
for some values of QmO using simulated data for future-SN 
survey like JDEM. Although we study a generic class of dark 
energy, e.g the thawing models, this shows the advantage of 
the large scale surveys like JDEM to distinguish different 
dark energy models from ACDM which is difficult with the 
presently available data. 

Next we study how effective is the simulated data 
in distinguishing different dark energy models using the 
Statefinder Hierarchies. In section 2, we already mentioned 
that S3 is the most promising parameter to distinguish dif- 
ferent dark energy models. Keeping this in mind, we calcu- 
late the bound on the present day value of 5*3, using the 
simulated JDEM data for different values of S7mo- In figure 
11, the regions inside the horizontal lines show this bound at 
2 — cr confidence level. This bound on S3 is valid for all thaw- 
ing models with power law potentials and is independent of 
the power law form. As one can see, the bound becomes 
narrower as one increases the QmO- Moreover, it is clear 
that as one decreases the value of SlmO, e.g. QmO = 0.24, 
53(2 = 0) = 1 is outside the 2 — cr allowed region ruling out 
the ACDM. But for higher values of flmo, ACDM is allowed. 
This is consisnent with our earlier conclusions from Fig 7 - 
Fig 10. 

In the same figure, we also show the behaviour of 



5*3(2 = 0) as a function of Ai for different potentials. For 
f2mO = 0.24, although ACDM is ruled out, but all the thaw- 
ing models more or less are inside the allowed region of S3 
making them difficult to distinguish. As one increases the 
QmO , models with higher A; start going out of the allowed re- 
gion of S3. This helps to distinguish different thawing mod- 
els with higher A^. We should emphasize the fact that for 
smaller Ai, all the thawing models have exactly similar be- 
haviour which is same as ACDM. It is only for higher Ai, 
the different thawing models behaves differently. 



5 DISCUSSION 

We look at observational constraints on the thawing scalar 
field models which have been proposed to explain late time 
cosmic acceleration. We first study the degeneracies of dif- 
ferent thawing models with other popular dark energy mod- 
els using the recently proposed ' Statefinder Hierarchy Sn 
'. We show that as one starts probing the higher order S„, 
different thawing models not only become distinguishable 
among themselves, they also behave significantly different 
from other dark energy models including ACDM. With ex- 
pected future data at significantly higher redshifts, this can 
be a smoking gun to distinguish different dark energy mod- 
els. We also study the behaviour of different thawing mod- 
els in the w' — w phase plane and compare those with the 
previously obtained LC bound. Using various observational 
data related to background evolution of the universe, we 
constrain the deviation of the thawing model from ACDM. 
With current precision level of the observational data ( at 
least for the background evolution), the models cannot be 
distinguished from the ACDM. But with simulated data for 
large future surveys like JDEM, we show that for some val- 
ues of f2mO, the data may distinguish thawing models with 
ACDM. Moreover, with future data, the S3 parameter may 
be quite useful in distinguishing different thawing models 
among themselves for certain plausible values of Q.mO- 

As a possible extension of this work, it will be really 
interesting to do a similar study with freezing type of scalar 
fields. As we mention in the introduction, Dutta and Scher- 
rer ( Dutta & Scherrer (2011)) have recently studied the 
slow-roll freezing scalar field models. In such a setup, one 



8 



can calculate the 'Statefinder Hierarchy' and compare them 
with those for the thawing models. With future high preci- 
sion data, such study can be useful to check which kind of 
scalar field model is actually preferred by the data. This will 
be our future goal. 
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